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We report on the realization of a stable solid state room temperature source for single photons. It 
is based on the fluorescence of a single nitrogen-vacancy (NV) color center in a diamond nanocrys- 
tal. Antibunching has been observed in the fluorescence light under both continuous and pulsed 
excitation. Our source delivers 2 x 10* s~^ single-photon pulses at an excitation repetition rate 
of 10 MHz. The number of two-photon pulses is reduced by a factor of five compared to strongly 
attenuated coherent sources. 

PACS numbers: 42.50.Dv, 03.67.-a, 78.67.-n 



I. INTRODUCTION 

The security of quantum cryptography is based on the 
fact that quantum mechanics does not allow one to du- 
plicate an unknown state of a single quantum system (for 
a review see This property is referred to as the no- 
cloning theorem. After the pioneering experiment of the 
group of Bennett and Brassard [|| , several quantum key 
distribution set-ups using attenuated laser pulses have 
been demonstrated (see for example refs [Q, §). In these 
implementations, single photons are approximated by 
strongly attenuated coherent pulses so that the average 
number of photons per pulse is pi « 0.1. In this case, 
the probability of having two photons in a pulse is ap- 
proximately Two-photon pulses are a potential 
information leakage source Indeed, an eavesdropper 
could tap on the communication between the sender (Al- 
ice) and the receiver (Bob), by keeping one of the two 
photons and sending the other one to Bob. The eaves- 
dropper can then measure the state of the photon once 
Alice and Bob have revealed their measurement basis. 
With attenuated coherent pulses, the only way to reduce 
the probability of having two photons in a pulse is to 
lower pi and thereby decreasing the transmission rate. 
Our aim is to realize an efficient single-photon source 
that would have a vanishing two photon probability for 
a non vanishing transmission rate. 

Single photons on demand can be produced by pulsed 
excitation of a single dipole ||, 0. The principle is that 
a single emitting dipole has to undergo a full excitation- 
emission-reexcitation cycle before emitting a second pho- 
ton. For a sufhciently short and intense excitation pulse, 
a single dipole emits one and only one photon [Q. 

Several solid state sources, like single organic molecules 

|, [10, 111, |l3| . self- assembled semiconductor quan- 
tum dots ||l4" Tsf " ^ , or semiconductor nanocrystals ||l7|] 
have been presented lately as potential candidates. How- 
ever, the ideal candidate should be photostable, work at 
room temperature and be easy to manipulate. Single 
nitrogen- vacancy (NV) color centers in bulk diamond |Tq] 



have recently been shown to exhibit strong antibunching 
at room temperature jl^ , [2C| , plf . They are intrinsically 
photostable and are believed to have a unity quantum ef- 
ficiency |2^. The high refractive index of diamond leads 
however to small a collection efficiency owing to total 
internal reflection and spherical aberrations. Also, the 
signal to background ratio is limited by the light emit- 
ted from the surrounding diamond crystal. One can also 
point out that bulk diamond cannot be manufactured in 
any desirable shape, and thus is very difficult to handle 
and insert, for example, into a microcavity. 

In this paper, we show that the use of single NV color 
center in diamond nanocrystals (typical size 50 nm) is 
a very convenient solution to these problems. The sub- 
wavelength size of the nanocrystals renders refraction ir- 
relevant. One can simply think of the nanocrystal as a 
point source emitting light in air. Furthermore, the small 
volume of diamond excited by the pump light reduces the 
emitted background light. Also, diamond nanocrystals 
can be easily handled in order to be inserted in a cavity 
or to be deposited on a fiber tip |^ . In addition NV cen- 
ters in diamond nanocrystals preserve all the important 
features of NV centers in bulk diamond. In particular, 
we have checked that the fiuorescence spectrum of NV 
centers in nanocrystals at room temperature is the same 
as in bulk. By investigating the autocorrelation function 
under continuous wave (CW) and pulsed excitation, we 
demonstrate the possibility to use NV centers in diamond 
nanocrystals as a room temperature stable single-photon 
source. 



II. EXPERIMENTAL SET-UP 

The color center used in our experiments is the NV 
defect center in synthetic lb diamond, with a zero phonon 
line at a wavelength of 637 nm jl^. The defect consists 
in a substitutional nitrogen and a vacancy in an adjacent 
site. A simplified level structure is a four-level scheme 
with fast non radiative decays within the two upper states 
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FIG. 1: Experimental set-up. The sample is excited by ei- 
ther a continuous or a pulsed frequency doubled YAG. Its 
fluorescence is collected using a confocal microscope set-up. 
The intensity correlations are measured using two avalanche 
photodiodes on each side of a 50/50 beam splitter, a time to 
amplitude converter and a multichannel analyzer. 
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and within the two lower states. This amounts to an 
incoherent two-level system. The lifetime of the excited 
state in the bulk is T{, = 11.6 ns 

NV centers are artificially created in synthetic MDA 
diamond powder from de Beers. Nitrogen is naturally 
present in diamond. Vacancies are created by irradia- 
tion with 1.5 MeV electrons at a dose of 3 x 10^^e~/cm^. 
Subsequent annealing in vacuum at 850°C during 2 hours 
leads to the formation of NV centers |l^ . The nanocrys- 
tals are dispersed by sonification in a solution of polymer 
(Polyvinylpyrrolidone at 1 weight% in propanol). This 
allows the disaggregation of the particles and their sta- 
bilization in a colloidal state. Centrifugation at 11000 
rpm for 30 min allows us to select nanocrystal sizes of 
do = 90 ± 30 nm (measured by dynamical light scatter- 
ing) . The average number of NV centers in a nanocrystal 
has been evaluated to 8 (see below). Nanocrystals con- 
taining a single NV center should therefore have a size 
around do/2 = 45 nm, which lies in the lower wing of the 
size distribution. The nanocrystal solution is then spin 
coated at 3000 rpm on thin fused silica substrates. Evap- 
oration of the solvent leaves a 30 nm thick film of poly- 
mer with the nanocrystals well dispersed on the surface. 
Their density was estimated to be around 0.02 fim^^. In 
most experiments we look at the centers from the other 
side of the plate, which is in contact with the oil of an 
immersion microscope lens (Nachet 004279, N.A. = 1.3). 

The experimental setup (fig. |l]), based on a confo- 
cal microscope, has been described in detail elsewhere 
^djj . The pump laser (CW or pulsed) is focused 
with an immersion oil, high numerical aperture objec- 
tive {NA = 1.3), onto a diffraction limited spot (« 400 
nm FWHM). The fluorescence light is collected by the 
same objective and after proper frequency and spatial 
filtering, it can either be send on a spectrometer, or to 
a Hanbury-Brown and Twiss set-up using two avalanche 
photodiodes (APD) from EG&G. Appropriate data pro- 



FIG. 2: (a) Confocal microscopy raster scan (5x5 fJ-ra^) of 
a diamond nanocrystal containing a single NV center. The 
count rate corresponds to one detector only. The size of a 
pixel is 100 nm and the integration time per pixel is 32 ms. 
The laser intensity impinging on the sample is 2.7 mW. In (b) 
a linescan along the dotted line of (a) is displayed, together 
with a gaussian fit, which is used to evaluate the signal and 
background levels. Here we obtain S/B — 20. Note that the 
fluorescence spot is slightly larger (500 nm) than the size of 
the excitation laser spot (400 nm) owing to saturation of the 
emitter. 

cessing allows us to obtain the histogram of the time sep- 
arations between successive photons. A slow (8s response 
time) x-y-z computerized servo-loop is used to compen- 
sate for any drifts. Bandpass filters allow detection from 
630 to 800 nm. This spectral window matches the broad 
emission spectrum of a NV center Q. 



III. CONTINUOUS WAVE EXCITATION 

At first, we investigate the fluorescence of single NV 
centers in nanocrystals under CW excitation with a fre- 
quency doubled YAG laser (A = 532 nm). Fig. |(a) 
displays a 2D scan of a single NV center. From the line 
scan in fig.^(b) we deduce a spatial resolution of 500 nm 
and a signal (S) to background (B) ratio of S/B — 20. 
The corresponding value in the bulk is = 8 (see fig. 

Fig. ^ shows the saturation behavior of the fluores- 
cence rate of NV centers in bulk and in nanocrystals with 
respect to pump power. The slightly decreasing rate at 
high pump power is due to a trapping metastable state 
[|20| . The saturating count rate in a diamond nanocrystal 
(•^nc ^ 4 4 X iQi g-i) is slightly lower than that in bulk 
diamond {N^ = 6.4 x 10'' s~'), but one has to take into 
account the longer lifetime of a NV center in a nanocrys- 
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FIG. 3: Fluorescence rate of a NV center in a nanocrystal (a) 
and in bulk diamond (b) as a function of the pump power. 
The count rate corresponds to one detector only. The crosses, 
empty squares, and black circles represent the background B, 
the total count rate T — S + B, and the signal from the 
center S — T — B, respectively. The data for the nanocrystal 
corresponds to the same center as that shown in fig. 0. 



tal (see at the end of this section). The number of pho- 



tons detected in a lifetime is TnrNV-'^ = 11 x 10 



the 



nanocrystal (Ufetime t„c = 25 ns) and TbN^ = 7.4 x 10^^ 
in bulk diamond (lifetime Tf, = 11.6 ns). This means that 
the geometrical collection efficiency for nanocrystals is 
increased by 50%. 

It can also be seen in fig. ^ that the contribution of 
background next to the nanocrystal is greatly reduced. 
The background coming from the diamond nanocrystal 
itself is reduced mainly because the excited volume of 
diamond is smaller. 

The evaluated overall detection efficiency is given by 
= vleoVoptVdet, whcrc i — b,nc for the bulk or for a 
nanocrystal, respectively. The geometrical collection ef- 
ficiencies for a dipole oriented orthogonally to the optical 
axis are rjg^^ = 0.18 and rjg^g = 0.38 (calculated). The 
optical transmission from the sample to the detectors is 
Vopt = 0.25 (measured), and the detector quantum ef- 
ficiency is rjdet — 0-7 (taken from the data sheet). In 
addition the NV center presents a metastablc state in 
which the excitation can be shelved [|l8|, |2^, |l]. This 
leads to a reduction of the count rate by a factor equal 
to the saturated population crf° — 0.25 ± 0.05 of the ex- 
cited state. This value has been inferred by fitting a 
saturation curve (cf fig. ||) and a set of autocorrelation 
functions (like fig. ^ obtained for different pump powers 
po| . This fit involves many parameters and gives only 
approximative results. For a nanocrystal, the saturated 
count rate should then be S"^ = ri^^a^/rnc = 6.6 x 10^ 
s"^, which is 15 times more than what we actually de- 
tected {Ns = 4.4 X lO'' s~^). This discrepancy exists also 
in bulk and its origin is still under investigation p2|. 
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FIG. 4: Autocorrelation function g^^^r) (left) and raw coinci- 
dence rate (right) for a NV center in a nanocrystal. The time 
bin w = 0.3 ns, total integration time is 323 s and the laser 
intensity impinging on the sample is 2.7 mW. Count rates on 
each photodiode are A^i = 22500 s"^ and N2 = 24500 s"^ 
The actual number of coincidences is indicated on the right. 
The zero-time value of the uncorrected normalized correla- 
tion function given by eq. (0) is Cn(0) = 0.17. The fit is 
performed with the model used in ]20(] . The data corresponds 
to the same center as that shown in fig. ^ 



We have also studied the autocorrelation function 
of the fluorescence of single NV centers in diamond 
nanocrystals. The raw coincidences c(r) (right axis) and 
autocorrelation function g'^^^r) = (/(t)/(i + T))/{I{t)Y 
(left axis) are represented in fig. |^. 

For evaluating the intensity correlation function g^^^ (t) 
of the NV center, the raw correlation data c(t) is normal- 
ized and corrected in the following way. The raw coin- 
cidence rate c(r) counted during a time T within a time 
bin of width w is first normalized to that of a Poissonian 
source according to the formula 



Cr(T) 



c{t)/{NiN2wT) , 



(1) 



where 7Vi_2 are the count rates on each detector. The 
normalized coincidence rate C™(t) is then corrected for 
the background light _B, and we obtain 
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(2 



)(r) = [C^-(r)-(l-p2)]/p2 



(2) 



where p = S/{S + B) is related to the signal to back- 
ground ratio, which is measured independently in each 
experimental run by measuring the count rate next to 
the nanocrystal (see fig. |^(b)). This takes into account 
only the background coming from the fused silica and 
the polymer, and not the parasitic light emitted by the 
diamond nanocrystal itself. Note that we have checked 
experimentally that the background light has a Poisso- 
nian statistics. 

It can be seen in fig. I that 5(2^(0) = 0.13, where 
the slight difference with zero is attributed to remaining 
background light emitted by the nanocrystal and elec- 
tronic jitter of the avalanche photodiodes (300 ps). This 
almost vanishing value of g^^^ (0) is the signature of the 
presence of a single emitter in the observed nanocrys- 
tal. In the case of the presence of p centers within a 
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nanocrystal, the value of the zero-time antibunching is 
1 — 1/p. This is actually how we estimate the number of 
NV centers in a nanocrystal. 



We have obtained a t — normalized coincidence rate 
(0) = 0.17 at the fluorescence rate maximum (input 
power of 2.7 mW), where the best value in bulk diamond 
was 0.26 As we shall see in the next section, 

this uncorrected normalized coincidence rate C™(0) is 
the relevant parameter for characterizing a single photon 
source. 

It should also be mentioned that g'^^^ (r) reaches values 
greater than unity for | r |> 10 ns. This bunching effect 
for longer time scale is due to the presence of the trap- 
ping metastable state in which the system can be shelved 
0, 1^, H Hi. This effect can also be seen as blinking 
on a time scale of « 50 ns. The time distribution of 
photons can ve viewed as bursts of photons of about 50 
ns duration. In each burst the delay between successive 
photons is always larger than 5 ns (antibunching). 

In a low pump regime, the central dip in the antibunch- 
ing traces can be fitted by an exponential function with 
the argument —T | r |, where F = j + r, with 7 being 
the NV center spontaneous decay rate and r the pump- 
ing rate [Q. Such fits have been performed for different 
pumping powers. The inverse lifetime 7 = 1/t„c of an 
NV center in a nanocrystal can then be inferred by ex- 
trapolating the value of the time constant for vanishing 
pump power (fig. We deduce a lifetime for NV cen- 
ters in bulk diamond of Tf, = 11.6 ± 0.1 ns ||2^ whereas 
the lifetime is found to be t„c = 25 ± 4 ns in diamond 
nanocrystals ||2^. This value has been obtained by ob- 
serving 10 different nanocrystals. A possible explanation 
for this lifetime increase is that the refractive index expe- 
rienced by the emitted light is different in bulk diamond 
and nanocrystals. When working out the spontaneous 
emission rate from the Fermi's golden rule, it turns out 
that this rate is proportional to the refractive index n of 
the material in which the dipole is radiating. 

In our case, the NV center in bulk diamond emits 
within a medium of index Ud = 2.4, whereas the cen- 
ter in a sub-wavelength nanocrystals emits as if it were 
in air for one half of the space, and in fused silica (refrac- 
tive index Ug — 1.45) for the other half. The expected 
lifetime is then r„c — Tf,[2nd/(1 -I- n^)] = 22.7 ns in good 
agreement with the experimental values. A full descrip- 
tion of lifetime changes due to refractive index is a con- 
troversial subject mainly because of local field correction 
issues ^,]3^ . However, our results tend to show that 
the local field experienced by the NV center in bulk and 
in nanocrystals is the same 27 1. As mentioned earlier 
no substantial difference in the emission spectrum of NV 
centers in nanocrystals and in bulk has been found. This 
is a good indication that the observed lifetime change is 
mainly due to the modification of the refractive index of 
the medium in which the NV center is radiating. 
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FIG. 5: Argument F of the exponential function fitting the 
antibunching traces obtained at different pump power. The 
extrapolation for vanishing pump power gives the lifetime. 
The black circles correspond to the data of fig. ^ and yield 
Tnc = 25 ns. The empty squares correspond to a NV center in 
bulk diamond (rj, — 11.6 ns). The slope for the nanocrystal 
is twice as large as that in bulk which is consistent with the 
lifetime increase, since the NV center in a nanocrystal can 
absorb twice as many photons during its lifetime. 



IV. PULSED EXCITATION 

The pulsed excitation consists of a home build source 
along the lines of reference |^ . The output of a 100 mW 
continuous wave NdiYAG laser (A = 1064 nm) is coupled 
into a fast (3 GHz) integrated LiNbOa modulator (Ale- 
nia) which slices up pulses of 1 ns duration at a repetition 
rate of 10 MHz. The pulses are then amplified to 1 W 
mean power by an Ytterbium fiber amplifier (Keopsys) 
and frequency doubled using a PPKTP crystal. In this 
way we obtain pulses at A = 532 ns of energy 2.5 nJ. 

We investigate the intensity autocorrelation function 
g^'^^r) of the fluorescence light of a single NV center in 
a nanocrystal under pulsed excitation. For a Poissonian 
source the probability of having a coincidence between 
two photons in the same pulse, or two photons coming 
from different pulses, is equal. Therefore the autocor- 
relation function for a pulsed Poissonian source exhibits 
peaks of same height separated by the repetition period 
(fig. i- 

In fig. ^ is shown the intensity autocorrelation func- 
tion of a single NV center under pulsed excitation. The 
excitation pulse duration is d = 1.2 ns and the repeti- 
tion period is 9 = 100 ns. It can be seen that the peak 
around r = is missing, which implies that the prob- 
ability of having two photons in one pulse is strongly 
reduced. This gives rise to a highly sub-Poissonian light 
source. In order to compare our single photon source to a 
pulsed Poissonian light source, one has to normalize the 
area c{m) of peak number m to the area N1N2OT of a 
peak corresponding to a pulsed Poissonian distribution 
with the same count rate, with T being the total acqui- 
sition time. Analogously to eq.(|l|), the normalized area 
of each peak is given by 

CAr(m) = c{m)/NiN2eT . (3) 

The determination of the area c(m) is performed in the 
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FIG. 6: Autocorrelation function from a fluorescing material 
(piece of white paper) excited by the pulsed laser. The emit- 
ted light has a Poissonian distribution. Pulse repetition pe- 
riod is 200 ns and the pulse width 1.2 ns. The count rates are 
Ni = 5011 and N2 = 5343 s"^ Integration time is T = 595 
s and the time bin is 2 ns. The area of each individual peak 
normalizes to unity. 
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FIG. 7: Pulsed autocorrelation function of a single NV cen- 
ter. Pulse repetition period is 100 ns, pulse width 1.2 ns and 
excitation mean power 0.9 mW. Count rates are A^'i = 10504 
s~^ and N2 = 9995 s~^. Integration time is T = 588 s and 
the time bin is 2 ns. The coincidences between peaks do not 
go down to zero because of the overlapping of adjacent peaks. 
The number above each peak represents its normalized area. 



following way. The peaks are fitted by exponential decays 
with the same lifetime. We checked that the lifetime 
found from the fit is the same than that deduced from the 
CW excitation. This fitting procedure allows an accurate 
evaluation of the area c(m) of peak to, in spite of the 
significant overlap between peaks. The normalized peak 
areas Cat (to) are given by the numbers displayed above 
each peak in fig. ^. For a Poissonian pulsed light source 
Cat (to) = 1 for all to (see fig. ||). For our single dipole 
in the most favorable case, in which the NV center is 
saturated but the zero time peak is as low as possible, we 
obtained Cn{0) = 0.21. Note that this value is slightly 
larger than C™(0) = 0.17. This is attributed to the finite 
duration {d = 1.2 ns) of the exciting pulses 0. 

Let us recall that the probability p2 of having two pho- 
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FIG. 8: Normalized peak area of fig. ^ as function of the peak 
number m. The experimental data is fitted with a function 
that assumes a random blinking of the NV center due to the 
trapping state. 



tons in a pulse is given by (assuming P2 ^ 1) 

P2 = CNiO)pi/2 (4) 

where pi is the probability of having a single photon. 
Note that for Poissonian light Cjv(O) = 1, and eq.(|4|) with 
Cn{0) = 1 gives the photon probability distribution of 
an attenuated coherent pulse. The zero time normalized 
coincidence rate Cn{0) = 0.21 means that the rate of two 
photon pulses is nearly five times lower than for Poisso- 
nian light. Since our source has a single photon rate of 
2 X 10'' s^^ at an excitation repetition rate of 10 MHz, 
its rate of two-photon pulses is only of 4 s~^. 

It can also be observed that the peaks m 7^ grow 
above unity. Just as in fig. U, this bunching effect for non 
zero time scale comes from the existence of a metastable 
state 0, m, |0[ |l[ |§ . In fig. | we have plotted the nor- 
malized area of the first peaks. The long time decay be- 
havior can be described by a simple model assuming that 
the NV center gets randomly trapped in the metastable 
state. This results in random blinking of the NV center. 
The normalized area of peak to is then given by [14| 



1 , Zk£p-(l/To// + l/To„)|m|( 
on 



(5) 



where To/ / is the mean time during which the excitation 
is trapped in the metastable state and the emission is in- 
hibited and Ton is the mean time during which the center 
is emitting. Fitting the normalized peak area (cf fig. |^) 
with eq.(^ allows us to extract the values Ton = 460 ns, 
Toff = 390 ns. The saturated count rate should then be 

7]^^[Ton/iTon + Toff)]/e = 1.4 X 10^ S'^, whilc 



^off 

ST 



we detect only iV^ = 2 x 10"* s~^. This is the same dis- 
crepancy that what was observed for CW excitation . 
Note that the factor Ton / {Ton+Tof f) = 0.54 accounts for 
the shelving state and plays the same role as = 0.25 
for CW excitation. Their different values are attributed 
to the fact that the shelving and deshelving rates depend 
on the excitation power pO| and are therefore different 
in CW and in pulsed regime. 
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V. CONCLUSION 

In this paper we have demonstrated the possibihty of 
using single NV centers in diamond nanocrystals as a 
room temperature sohd state source for single photons. 
They present the advantages of being photostable and 
easy to manipulate. Furthermore the fabrication of the 
samples is easy and inexpensive. The single photon rate 
summed over both photodiodes is 2 x 10^ s^^ for an ex- 
citation repetition rate of 10 MHz. Our source exhibits 
a strong sub-Poissonian distribution. The two photon 
pulse rate is reduced by a factor of five compared to a 
Poissonian source and is equal to 4 s~^. Improvements 
in the collection efficiency should be obtained by deposit- 
ing the nanocrystals on a mirror or inserting them into a 
microcavity. Even though the repetition rate is low com- 
pared to what can be obtained with attenuated pulses 



using laser diodes, the reduction by a factor of 5 of the 
two photon pulses will allow a secure transmission over 
larger distances 
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